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We present a systematic analysis at the leading log order of the influence of combination of bilinear 
and trilinear i?-parity violating couplings on the decay b — ► s + 7 . Such contributions have never 
been explored in the context of b — ► s + 7 decay. We show that influence of charged-slepton-Higgs 
mixing mediated loops can dominate the SM and MSSM contributions and hence can provide strong 
bounds on the combination of bilinear-trilinear i?-parity violating couplings. Such contributions are 
also enhanced by large tan /3. With substantially extended basis of operators (28 operators) , we 
provide illustrative analytical formulae of the major contributions to complement our complete 
numerical results which demonstrate the importance of QCD running effects. 

PACS numbers: 



Introduction. — The minimal supersymmetric standard 
model (MSSM) has been the most popular candidate the- 
ory for physics beyond the Standard Model (SM) for the 
last couple of decades. With the recent accumulation of 
evidence for neutrino oscillations, it is clear that the lep- 
ton number conserving MSSM has to be amended. The 
simplest option is then to give up imposing i?-parity and 
hence admit all gauge invariant terms in a (generic) su- 
persymmetric SM. Other alternatives include incorporat- 
ing any particular neutrino sector model such as the see- 
saw mechanism with extra gauge singlet superfields. We 
focus on the first optionQ, The model has the spe- 
cial merit that the parameters that give rise to neutrino 
masses and mixings also have interesting phenomenolog- 
ical consequence in the quark and charged lepton sectors. 
Here in this letter, we report on some of the novel con- 
tributions to the b — ► s + 7 decay from the model. 

Within the SM, the flavor sector still needs more 
scrutiny from theory as well as experiments. In partic- 
ular, flavor changing neutral current (FCNC) processes 
are widely considered to be the window of the physics 
beyond SM. Among the processes, B — > A" s + 7 is a par- 
ticularly attractive candidate. The most up-to-date SM 
prediction Q gives 

Br[B -^X s + 1{El > 1.6 oev)] SM =(3.57 ± 0.30) x 10~ 4 , 

. n (1) 

while the experimental number (world-average) is Q 

Br[B -> X s + 7 ( e 7 > 1.6 gcv)] exp =(3.34 ± 0.38) x 10~ 4 . 

.(2) 

It clearly leaves not much room for new physics contribu- 
tions. Hence, can be used to obtain stringent constraints 
on flavor parameters of various new physics models. 

There have been some studies on the process within the 
general framework of i?-parity violation 0, [(j. Ref.[j|, 
fails to consider the additional 18 four-quark operators 
which, in fact, give the dominant contribution in most 
of the cases. The more recent work of ref.[6| has consid- 
ered a complete operator basis. However, we find their 
formula for Wilson coefficient incomplete. In fact, the 
particular type of contributions — namely, one from a 
combination of a bilinear and a trilinear i?-parity violat- 



ing (RPV) parameters, we focus on here, has not been 
studied in any detail before. 

It is not possible for us to give much analytical details 
of our study here in this short letter. We will only outline 
the major features of the full analysis given in a paral- 
lel report @, to which interested readers are referred. 
We adopt an optimal phenomenological parametrization 
of the full model Lagrangian, dubbed the single single- 
VEV parametrization (SVP), first explicitly advocated 
in refs.|8|]. It is essentially about choosing a basis for 
Higgs and lepton superfields in which all the "sneutrino" 
VEVs vanish. The formulation gives the simplest expres- 
sions for all the mass matrices of the fermions and scalars 
without a priori assumption on the admissible form of im- 
parity violation. In particular, all the RPV effects to the 
fermion mass matrices are characterized by the three bi- 
linear parameter /Vs. Working under the formulation, 
it has been pointed out in ref.Q that there are interest- 
ing contributions to the down squark and charged slep- 
ton mass matrices of the form (A«*A' ijfc ) and (fj,*\ijk )■ 
These give explicit indications of the existence of bilinear- 
trilinear type contributions to fermion dipole moments at 
1-loop order, including the transitional moment term to 
be identified with the b — > s + 7 decay. Detailed analyt- 
ical and numerical studies have been performed on the 
case of neutron electric dipole moment (lfj and /i — > e + 7 
[ilf. Here, we report on the more difficult calculation 
of B — > X s + 7 . Background details on the model and 
the various mass matrices are given in ref.0. Based on 
the latter results, we implement our (1-loop) calculations 
using mass eigenstate expressions jjj, hence free from the 
commonly adopted mass-insertion approximation. While 
a trilinear RPV parameter gives a coupling, a bilinear pa- 
rameter now contributes only through mass mixing ma- 
trix elements characterizing the effective couplings of the 
mass eigenstate running inside the loop. The /Vs are 
involved in fermion, as well as scalar mixing sQ. There 
are also the corresponding soft bilinear Bi parameters 
involved only in scalar mixings^]. Combinations of fi^s 
and Bi's with the trilinear A' ijfc parameters are our major 
focus. 

The Effective Hamiltonian Approach. — The partonic 
transition b — > s + 7 is described by the magnetic pen- 
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guin diagram. Under the effective Hamiltonian approach, 
the corresponding Wilson coefficients of the standard Q 7 
operator has many RPV contributions at the scale M w . 
For example, we separate the contributions from differ- 
ent type of diagrams as C 7 = C) v + Cf + Cj + C£ + 

C-j + C7 . The neutral scalar loop contribution , 
as an illustrative case, is proportional to 



Qd 



F 3 



ML 



M L 



where the effective vertex couplings Af^ mi and Af^ mi for 
the mass eigenstates each contains the a A'-coupling 
contribution!?! ■ 

Apart from the 8 SM operators with additional con- 
tributions, we actually have to consider many more op- 
erators with admissible nonzero Wilson coefficients at 
M w resulting from the RPV couplings. These are the 
chirality-flip counterparts Q7 and Qs of the standard 
(chromo) magnetic penguins Q7 and Qg, and a whole list 
of 18 new relevant four-quark operator of current-current 
type to be given as : 



23,4 = (sRa l^bRa.p) ^ (QRiP 7^ 9ffi/3 
i— u.c.d.s.b 



d, s, 6; (3) 



25,6 = (s Ra ^b Ra! 0) 



i— u.c.d.s.b 



; (4) 

(q~Li(3i^ qhip.a) ; (5) 



2g-i3 = (sRa Y l b R p) {q L p 7m > q = d,s,b, u, c(6) 

and six more operators from A" couplings^ we skip here 
for brevity. The interplay among the full set of 28 oper- 
ators is what makes the analysis complicated. The effect 
of the QCD corrections proved to be very significant even 
for the RPV parts. 

We skip here the details involved in the evaluation of 
the various effective Wilson coefficients for the decay rate 
of b — > s + 7 and give only the numerical results from our 
leading log (LL) order analysis |7j : 



Cf(m b ) 



-0.351 Cf{M w ) + 0.665 Cf(M n 
-0.093 CJ ff (M w ) - 0.198 Cf(M n 
-0.198 C^(M W ) - 0.178 Cff(M w 



Cf(m b ) = 0.381 Cf(M w ) 
+0.093 Cf(M n 
-0.198 C^f(M w 
+0.510 Cff(M w 
+0.381 Clf(M n 



0.665 Cf(M w ) 

- 0.198 Cf{M w ) 

- 0.178 Ctf{M w ) 
+ 0.510 Glf{M w ) 

- 0.213 Ctf{M w ) . (7) 



The branching fraction for Br(b — 
through the semi-leptonic decay b - 
large bottom mass dependence (~ ra\ 
in CKM elements cancel out. 



s + 7) is expressed 
► u\cev so that the 
') and uncertainties 



Br(b - s + 7 ) = ffj^j*^) Br^ p (b - u\ce D e ) 



where Br exp (b — ► u\cev e ) — 10.5% and 



64tH 



\Cf(» b )\ 2 + \Cf{n b )t 



(8) 



(9) 



Note that we have also to include RPV contributions to 
the semi-leptonic rate for consistency .7]- 

Analytical Appraisal of the Results. — There are three 
kinds of bilinear RPV parameters, Hi,Bi and m| . re- 
lated by the tadpole equation constraints 0,|9j]. Without 
loss of generality, we choose and Bi to be independent. 
The influence of a \Bi\ (or in conjunction with 
is felt through the lepton number violating mass mixings 
in (s)leptonic propagators of tree and penguin diagrams. 
\Bi\ insertions may have much stronger influence than 
the \/j,i\ as the former case is inversely proportional to 
light slepton mass-squared whereas the latter ones come 
with the inverse of heavier squark mass-squared. We fo- 
cus our discussion here on the \Bi\ insertions to provide 
an analytical appraisal of the numerical results. The case 
for fii can be appreciated in a similar fashion [l3|. 

There are two kinds of B^-X 1 combinations that con- 
tribute to 6 — > s + 7 at 1-loop: (a) B*\^ 2 , and (b) 

Bi\* s . These involve quark-scalar loop diagrams. Case 
(a) leads to the b L —> s R transition (where SM and MSSM 
contribution is extremely suppressed) whereas case (b) 
leads to SM-like b R — > s L transition. For the purpose of 
illustration, we will assume a degenerate slepton spec- 
trum and take the sleptonic index i — 3 as a represen- 
tative. The j = 3 [2] contributions for case (a) [(b)] 
with both sneutrino-Higgs mixings and charged-slepton- 
Higgs mixings are easy to appreciate. For the j values, 
the charged loop contributions are still possible by invok- 
ing CKM mixings. Consider the contribution of case (a) 

with I-B3A332I to the Wilson coefficient Cj, for instance. 
Through the extraction of the bilinear mass mixing ef- 
fect under a perturbative diagonalization of the mass 
matrices we obtain 



Cf 



M? 



y b tan/3 
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(10) 
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for the charged and neutral scalar loop, respectively. In 
the above equations, proportionality to tan/3 shows the 
importance of these contributions in the large tan/3 limit. 
The Mg, M/, M s 2 , are all scalar (slepton/Higgs) mass 
parameters. The term proportional to y t above has chi- 
rality flip into the loop. Thinking in terms of the elec- 
troweak states, it is easy to appreciate that the loop dia- 

gram giving a corresponding term for (cf. involving 
A/J m3 A^ 2 ) requires a Majorana-like scalar mass inser- 
tion, which has to arrive from other RPV couplings 0. 
In the limit of perfect mass degeneracy between the scalar 
and pseudoscalar part (with no mixing) of multiplet, it 
vanishes. Dropping this smaller contribution, together 
with the difference among the Inami-Lim loop functions 
and the fact that the charged loop has more places to at- 
tach the photon (with also larger charge values) adding 

up, we expect the Cj to be larger than C$ . 

Numerical Results. — The basic strategy for the im- 
plementation of our numerical study is to choose model 
parameters in such a way that in the limit that the 
RPV parameters vanish, the resulting MSSM gives the 
decay rate well within the experimental limit 12]. The 
otherwise arbitrariness in the choice mostly does not 
have much of an effect on the qualitative dependence of 
the results on a specific combination of RPV parame- 
ters. Note that we always keep i?-parity conserving fla- 
vor violating squark and slepton mixings vanishing, to 
focus on the RPV effects. We take non-vanishing val- 
ues for relevant combinations of a bilinear and a trilin- 
ear RPV parameters one at a time, and stick to real 
values only. Our model choice is (with all mass dimen- 



sions given in GeV): squark masses 300, down- type Higgs 
mass 300, [i = —300 sleptons mass 150 and gaugino 
mass M 2 = 200 (with Mi = 0.5M 2 and M 3 = 3.5M 2 ), 
tan/3 = 37 and A parameter 300. The mass for H u and 
soft bilinear parameter Bq are determined by electroweak 
symmetry breaking conditions which are modified in the 
presence of RPV parameters [lj, || . 

Under the scenario discussed, we impose the experi- 
mental number to obtain bounds for each combination 
of RPV parameters independently (given in Table I) . We 
address here a couple of cases in a bit more detail. Con- 
sider, for instance, the case (b) combination f? 3 A 3 * 3 |. We 
obtain a bound of 5.0 x 10 -5 , when normalized by a factor 
of /J, 2 . Since this is a 6 H — ► s L transition, the RPV contri- 
bution interferes with the SM as well as the MSSM con- 
tribution. In Fig. n we have plotted the relevant Wilson 
coefficients. Over and above the loop contributions we 
see that there are contributions coming from four-quark 
operator with Wilson coefficients C\\ (oc yb) which is 
stronger than the other two four-quark quark coefficients 
Cio,i3 oc y s (not shown in the graph). Since the neutral 
scalar loop contribution is proportional to the loop func- 
tion F\ (which is or order .01), it is suppressed compared 
to current current contributions. Also here the charged 
scalar contribution comes only with chirality flip inside 
the loop and has a CKM suppression. So the current- 
current is dominant. It has a more subtle role to play 
when one writes the regularization scheme-independent 
Cy ff = Cf — Cxi at scale M w . Due to dominant and neg- 
ative sign chargino contribution (A t /x < 0), the positive 
sign Cn interferes constructively with C7 and enhances 
the rate. 

The case (b) combination I-B3A332I is a different story, 
as it leads to b L — > s R transition and hence RPV does not 
interfere with SM or MSSM contribution. This leads to a 
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FIG. 1: Various Wilson coefficients versus |i?3 A 32 3 1 ■ '+' sign 
stands for MSSM chargino contribution, 'x' stands for the 
MSSM charged Higgs contribution, '★' stands for sneutrino 
contribution, all contributing to C7. 'Empty square' stands 
for Cn(Mw), 'filled square for C~i(Mw) and 'empty circle' 
for C7(m,b). 
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FIG. 2: Various Wilson coefficients versus 
stands for charged slepton contribution, 'x 
trino contribution, '★' stands for MSSM chargino contribu- 
tion, all contributing to C7. 'Empty square' stands for MSSM 
charged Higgs contribution, 'filled square' for Cu(Mw), 
'empty circle' for Cv(M w ), 'filled circle' for CV(Mw), 'empty 
triangle' for CV(mf,), and 'filled triangle' for CV(mb). 
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TABLE I: Bounds for the products of bilinear and tri- 
linear RPV couplings. 



Product 


Our bound 


Wilson Coef. 


g ' X j 2a 

Mil 




5.0 x 10 -5 


CV.S, 6*7,8, C\l, C10, C13 


S i' A i32 




7.4 x 10" 3 


C7,8, CV,8, C10, C'll 






2.3 x 10" 3 


CV,8, CV,8 


Mn 




6.5 x 10" 2 


C7,8 


»l 




8.0 x 10 -2 


CV,8 


S i' A il2 




4.5 x 10" 2 


C7,8 


Mi' A i23 
MO 




2.2 x 10 -3 


C7,8 


Mi' A i32 
MO 




1.0 x 10" 2 


CV,8 


Mi ' A i33 
MO 




8.0 x 10" 2 


C7,8 



bound of 7.3 x 10~ 3 after normalization by In Fig. [3 
we have again plotted the relevant Wilson coefficients. 
Unlike the previous case, here we see that the contri- 
butions from charged scalar loop dominates over both 
neutral scalar as well as current-current contributions. 
This is in accordance with our analytical expectations. 
Again, the current-current contribution due to G\i has a 
very subtle role to play here. The regularization scheme- 
independent effective Wilson coefficient Cf = C7 — C\\ 
at the scale M w . The negative sign leads to cancellations 
and hence weakens the bound. 

The influence of |MiA^ fe | is less stronger than the Bi in- 
sertions. The Hi insertions affect the MSSM chargino and 
the neutralino type of diagrams by mixing them with the 
charged and the neutral leptons. Since such contributions 
are suppressed by the heavier squark masses, the influ- 
ence is not as strong as slepton loops. However there exist 
gluino mediated loop diagrams with a flavour violating 
chirality flip in the down-squark propagator (cx |/J.iA^ fe |) 
which gives non-negligible contributions and indeed lead 
to good bounds. 

Conclusions. — To conclude we have systematically 



studied the influence of the combination of bilinear- 
trilincar RPV parameters on the decay b — > s + 7 an- 
alytically as well as numerically. Such a study has not 
been attempted before. We demonstrate the plausible 
dominance of the RPV contributions over conventional 
SM and MSSM contributions in some parameter space 
regions. These contributions are enhanced by large tan/3. 
It is shown that charged-slepton-Higgs mixing mediated 
loop typically dominates over the sneutrino-Higgs mix- 
ing mediated loop. Our study has consistently incorpo- 
rated all the QCD corrections at the leading log order, 
with a whole list of extra operators and their Wilson co- 
efficients arising from RPV couplings. We have shown 
that, through the formulation of scheme-independent ef- 
fective Wilson coefficients, the new current-current oper- 
ators can considerably influence the decay rate. Under a 
typical and compatible model parameter choice, we ob- 
tain strong bounds on several combinations of RPV pa- 
rameters. Bounds on such bilincar-trilinear parameter 
combinations are not available before. Various b — » s + 7 
contributions over different parameter space regions may 
complicate the story and partial cancellation among them 
are a likely possibility. And our leading log calculation 
bears relatively large certainty. Nevertheless, the bounds 
show values of the RPV parameter combinations that will 
play a major role in endangering the compatibility of the 
theoretical b — ► s + 7 result with the experimental limits. 
This interpretation of our results is quite robust. 

Our analytical formulae include all RPV contributions 
at 1-loop level. Numerical study has also been performed 
on combinations of trilinear parameters 0. We quote 
here a few exciting bounds under a similar sparticle spec- 
trum. For instance |A^ 33 • A^l for i = 2, 3 should be less 
than 1.6 x 10~ 3 to be compared with rescaled existing 
bound of 2 x 1CT 2 . 
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